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INITIAL-BOUNDARY VALUE PROBLEM FOR 
INTEGRABLE NONLINEAR EVOLUTION EQUATIONS 
WITH 3x3 LAX PAIRS ON THE INTERVAL 

JIAN XU* AND ENGUI FAN 

Abstract. We present an approach for analyzing initial-boundary 
value problems which is formulated on the finite interval (0 < a: < 

L, where L is a positive constant) for integrable equations whose 
Lax pairs involve 3x3 matrices. Boundary value problems for in¬ 
tegrable nonlinear evolution PDFs can be analyzed by the unified 
method introduced by Fokas and developed by him and his collab¬ 
orators. In this paper, we show that the solution can be expressed 
in terms of the solution of a 3 x 3 Riemann-Hilbert problem. The 
relevant jump matrices are explicitly given in terms of the three 
matrix-value spectral functions s{k),S{k) and S'i(fc), which in turn 
are defined in terms of the initial values, boundary values at x = 0 
and boundary values at x = L, respectively. However, these spec¬ 
tral functions are not independent, they satisfy a global relation. 

Here, we show that the characterization of the unknown boundary 
values in terms of the given initial and boundary data is explicitly 
described for a nonlinear evolution PDF defined on the interval. 

Also, we show that in the limit when the length of the interval 
tends to infity, the relevant formulas reduce to the analogous for¬ 
mulas obtained for the case of boundary value problems formulated 
on the half-line. 


1. Introduction 

Integrable PDEs have the distinctive property that they can be writ¬ 
ten as the compatibility condition of two linear eigenvalue equations, 

Date: September 10, 2015. 

Key words and phrases. Riemann-Hilbert problem, Manakov Systems, Initial¬ 
boundary value problem. 


1 



2 


J.XU AND E.FAN 


which are called a Lax pair [T]. An effective method (Inverse Scatter¬ 
ing Transform (1ST)) for solving the initial value problem for integrable 
evolution equations on the line was discovered in 1967 |2]. However, 
the presence of a boundary presents new challenges. It was realized in 
[3] that the extension of this method to initial boundary value problems 
requires a deeper understanding of the following question: What is the 
fundamental transform for solving initial boundary value problems for 
linear evolution equations with x-derivatives of arbitrary order? The 
investigation of this question has led to the discovery of a general ap¬ 
proach for solving boundary value problems for linear and for integrable 
nonlinear PDEs [1] (see also [HE]). For integrable nonlinear evolution 
PDEs this approach is based on the simultaneous spectral analysis of 
the two linear eigenvalue equations forming the Lax pair, and on the 
investigation of the so-called global relation, which is an algebraic re¬ 
lation coupling the relevant spectral functions. 

The Fokas method provides a generalization of the 1ST formalism 
from initial value to initial-boundary value (IBV) problems, and over 
the last eighteen years, this method has been used to analyze boundary 
value problems for several of the most important integrable equations 
with 2x2 Lax pairs, such as the Korteweg-de Vries |7j, the nonlinear 
Schrodinger [S], the sine-Gordon equations [9], see [iniEiDsiE]. Just 
like the 1ST on the line, the unified method yields an expression for 
the solution of an IBV problem in terms of the solution of a Riemann- 
Hilbert problem. In particular, the asymptotic behavior of the solution 
can be analyzed in an effective way by using this Riemann-Hilbert 
problem and by employing the nonlinear version of the steepest descent 
method introduced by Deift and Zhou HZI. 

In 2012, Lenells first develops a methodology for analyzing IBV prob¬ 
lems on the half-line for integrable evolution equations with Lax pairs 
involving 3x3 matrices [IB]. Although the transition from 2 x 2 to 
3x3 matrix Lax pairs involves a number of novelties, the two main 
steps of the method of mM remain the same: 
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(1) Construct an integral representation of the solution characterized 
via a matrix Riemann-Hilbert problem formulated on the complex k- 
sphere, where k denotes the spectral parameter of the Lax pair. This 
representation involves, in general, some unknown boundary values, 
thus the solution formula is not yet effective. 

( 2 ) Characterize the unknown boundary values by analyzing the so- 
called global relation. In general, the characterization of the unknown 
boundary values involves the solution of a nonlinear problem. 

After Lenells’ work, IBV problems on the half-line for other inte- 
grable evolution equations such as the Degasperis-Procesi [IH], Sasa- 
Satsuma [20] , three wave [ 21 ] , the two-component nonlinear Schrodinger 
[22] . the Ostrovsky-Vakhnenko [22] equations, are analyzed. However, 
within the knowledge of the authors, the IBV problems for integrable 
equations with 3x3 matrices Lax pair on the hnite interval has not 
been studied yet. 

The purpose of this paper is to extend the ideas from analyzing the 
IBV problems on the half-line to the hnite interval for integrable evoul- 
tion equations with Lax pairs involving 3x3 matrices. In fact, dealing 
with IBV problems on the interval has some difficulties. The imple¬ 
mentation of step ( 1 ), we need four curve integration from the four 
corners of the (x, f)—domain. We will dehne analytic eigenfunctions, 
denoted by {M„(x, f, A;)}, via integral equations which involve integra¬ 
tion from all four corners simultaneously. The most difficulties is to 
make a distinction between the integration contour 73 and 74 when we 
try to analyze the IBV problems on the interval. It is different from 
the analyzing the IBV problems on the half-line, because in the half- 
line case there just one integration curve 73 . Here, the constructions 
of this paper can be compared with the corresponding formalism for 
2 X 2-matrix Lax pairs introduced by Fokas and Its, see m- The im¬ 
plementation of step ( 2 ), the differences are introducing a new factor 

during analyzing the global relation to characterize the unknown 
boundary data in terms of the given initial and boundary data. We 
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show that in the limit when the length of the interval tends to inhty, 
the relevant formulas reduce to the analogous formulas obtained for 
the case of boundary value problems formulated on the half-line. 

Organization of the paper: In the next subsection, we introduce 
our main example considered in this paper. In section 2 we perform 
the spectral analysis of the associated Lax pair. We formulate the 
main Riemann-Hilbert problem in section 3 and this concludes the 
implementation of step (1) above. We also get the map between the 
Dirichlet and the Neumann boundary problem through analyzing the 
global relation in section 4 and this concludes the implementation of 
step (2). 


1.1. The main example. In this paper, we will consider the two- 
component nonlinear Schrodinger equation or Manokov equation 


iqit + Qixx - 2a(|gip -f 1^21^)91 = 0, 

iq2t + q2xx - 2a(|gip -f 1^21^)92 = 0. 


a = ± 1 . 


( 1 . 1 ) 


where qi{x,t) and q 2 {x,t) are complex-valued functions of (x,f) G fl, 
with denoting the hnite interval domain 


= {(a;,f)|0 < X < L,0 < t < T}, (1.2) 


here L > 0 is a positive hxed constant and T > 0 being a hxed hnal 
time. Here, a = 1 means defocusing case and a = —1 means focus¬ 
ing case. This system was hrst introduced by Manakov to describe 
the propagation of an optical pulse in a birefringent optical hber p4] . 
Subsequently, this system also arises in the context of multicomponent 
Bose-Einstein condensates [25] . 
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We will consider the following initial-boundary value problem for the 
2-NLS equation, 


Initial value: gio(T) = gi(T, f = 0), ? 2 o( 3 ^) = 92 (t, f = 0), 

Dirichlet boundary value: goi{t) = qi{x = 0, t), gmif) = Q 2 {x = 0, t), 

foiit) = qi{x = L,t), /o 2 (t) = q 2 {x = L,t), 

Neumann boundary value: gu{t) = qix{x = 0, t), gi 2 (t) = q 2 xix = 0, t), 

fn{t) = qix{x = L,t), /i 2 (t) = q 2 x{x = L,t). 

(1.3) 

It is well known that 2-NLS equation admits a 3 x 3 Lax pair, 




/ \ 

^2 

V ^3 / 


(1.4a) 


= V^. (1.4b) 

where 

U = ikA + Vi. (1.5) 

and 

V = 2iPA + V 2 (1.6) 


here 



/ -1 0 0 \ 


( ° 

Qi 0.2 \ 

A = 

) 0 

^ 0 

,h7 = 

aqi 

') 0 

) 0 


\ 0 0 1 / \aq2 0 0 / 


where 


2kV2^^^ + 

(1.7) 





= zA(V,^ - V,x). 


( 1 . 8 ) 


2. Spectral analysis 


2.1. The closed one-form. Introducing a new eigenfunction/^( t, f, A:) 
by 

^ = ^ei^kx+2iKkH (^ 2 . 1 ) 
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then we find the Lax pair equations 

f - [ikK^n] = Vin, 

[fit - [2ik^A, fi] = V2fi. 

Letting A denotes the operators which acts on a 3 x 3 matrix X by 
AX = [A, X] , then the equations in fl2.2p can be written in differential 
form as 

= W, (2.3) 

where W(x, t, k) is the closed one-form defined by 

W = + V2dt)fi. (2.4) 


2.2. The fij’s definition. We define four eigenfunctions {fijjf of fl 2 . 2 p 
by the Volterra integral equations 

fij(x,t,k) = I+ f j = l,2,3,4. (2.5) 

J'yj 

where Wj is given by fl2.4p with fi replaced with fij^ and the contours 
{ffili showed in Figure 1. The first, second and third column of 


ti ti ti t 



o ^ o ^ o ^ o 


7i 72 73 74 

Figure 1. The three contours 71 , 72,73 and 74 in the (x, t)—domain. 


the matrix equation 

[AjA- 

[Ajh- 

[AjA'- 


fl2.5p involves the exponentials 

^ik{x—x^)+AikP'{t—t') 

^—2ik{x—x')—Aik‘^{t—t') 

^—2ik{x—x')—4ik‘^{t—t') 


( 2 . 6 ) 
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And we have the following inequalities on the contours: 

71 : T — x' > 0, t — t' < 0, 

72 : X — x'> 0 , t — > 0 , , , 

73 : X — x' < 0, t — t' > 0, 

74 : X — x' < 0, t — t' < 0. 

So, these inequalities imply that the functions {Hj}i are bounded and 
analytic for /c G C such that k belongs to 


k'l'■ {^2,03,03), 

fi 2 ■ (-Di,-D4,-D4), 

fJ'S ■ {D3,D2,D2), 

^4 : (-D 4 , Di, Di). 


( 2 . 8 ) 


where {Dn}\ denote four open, pairwisely disjoint subsets of the com¬ 
plex /c—sphere showed in Figure 2. 



Figure 2. The sets Dn, n = 1,..., 4, which decompose 
the complex A:—plane. 
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We also notice that the sets {Dn\\ has the following properties: 

Di = {k E C|ReZi > Re/2 = Re/3, Rezi > Re^2 = Re^s}, 

D2 = {k E C|Re/i > Re/2 = Re/3, Rezi < Re^2 = Re2:3}, 

Dz = {k E C|Re/i < Re/2 = Re/3, Rezi > Rez2 = Re2:3}, 

-D4 = {/s G C|Re/i < Re/2 < Re/3, Rezi < Rez2 = Re-2;3}, 

where li{k) and Zi{k) are the diagonal entries of matrices ikA and 2 //c^A, 
respectively. 


2.3. The M„’s definition. For each n = 1,...,4, dehne a solution 
k) of fl2.2p by the following system of integral equations: 


{Mn)ij{x,t, k) 





t\k)) 




k E Dn, z,j = 1,2,3. 
(2.9) 


where Wn is given by fl2.4|) with p replaced with Mn, and the contours 


7*", n = 

= 1,- 

..,4, 

i,j = 1 

,2 

3 are dehned 

by 




7i 

if 

Reli{k) 

< 

Relj{k) 

and 

Rezi{k) 

> 

Rezj{k) 

77 , 

72 

if 

Reli{k) 

< 

Relj{k) 

and 

Rezi{k) 

< 

Rezj{k) 

R, = ^ 

73 

if 

Reli{k) 

> 

Relj{k) 

and 

ReZi{k) 

< 

ReZj{k) 


. 74 

if 

Reli{k) 

> 

Relj{k) 

and 

Rezi{k) 

> 

Rezj{k) 


( 2 . 10 ) 

Here, we make a distinction between the contours 73 and 74 as follows. 



ni<i<j< 3 (Ile^i(^) - Relj{k)){ReZi{k) - Rezj{k)) < 0, 
ni<i<,< 3 (Ile/i(A:) - Relj{k)){ReZi{k) - Rezj{k)) > 0. 

( 2 . 11 ) 

The rule chosen in the produce is if Im = In, xa may not equals n, we 
just choose the subscript is smaller one. 
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According to the definition of the 7 ”, we find that 



( 

74 

74 

74 



( 

73 

73 

73 

\ 

7 = 


72 

74 

74 




7i 

73 

73 



\ 

72 

74 

74 



\ 

7i 

73 

73 

/ 


/ 

73 

7i 

7i 

\ 


/ 

74 

72 

72 

\ 

7^ = 


73 

73 

73 


74 = 


74 

74 

74 



V 

73 

73 

73 



V 

74 

74 

74 

/ 


( 2 . 12 ) 


The following proposition ascertains that the M„’s defined in this 
way have the properties reqnired for the formulation of a Riemann- 
Hilbert problem. 


Proposition 2.1. For each n = the function Mn{x,t,k) is 

well-defined by equation \2.9t) for k G Dn and {x, t) E fl. For any fixed 
point {x,t), Mn is hounded and analytic as a function of k E Dn away 
from a possible discrete set of singularities {kj} at which the Fredholm 
determinant vanishes. Moreover, Mn admits a hounded and contious 
extension to Dn and 


Mn(x,t,k) =I-\- 0 (y), A: —)■ 00 , k E Dn- (2-13) 

k 

Proof. The bounedness and analyticity properties are established in 
appendix B in [18]. And substituting the expansion 


M = Mo + + • • • , 


k 


k^ 


k ^ 00. 


into the Lax pair fl2.2l) and comparing the terms of the same order of 
k yield the equation (12.131) . □ 


2.4. The jump matrices. We define matrix-value functions Sn{k), 
n = 1,..., 4, and 


Sn{k) = Mn{0,0,k), k E Dn, 


n = 1 ,..., 4 . 


( 2 . 14 ) 
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Let M denote the sectionally analytic function on the complex A:—sphere 
which equals for k G Dn- Then M satishes the jump conditions 

-Thi k G Dn Ll I-5 * * * ; ^5 Tl ^ Tfl^ 

(2.15) 

where the jump matrices k) are dehned by 

= (2.16) 


2.5. The adjugated eigenfunctions. We will also need the analytic- 
ity and boundedness properties of the minors of the matrices f, k)}\. 

We recall that the cofactor matrix of a 3 x 3 matrix X is dehned 
by 




( mii(W) -mi 2 (W) mi 3 (X) ^ 

-m 2 i(X) m 22 (X) -m23(X) , 

\ 77131 (X) -77732 (W) 77733 (X) / 


where mij{X) denote the (ijjth minor of X. 

It follows from fl2.2p that the adjugated eigenfunction satishes 
the Lax pair 


hx + [ikA,n^] = 

fif + [2ik‘^A, fi^] = -Vj' /7^. 


(2.17) 


where V'^ denote the transform of a matrix V. Thus, the eigenfunctions 
{/7^}^ are solutions of the integral equations 


Af 


{x,t,k) = I - / + j = 1,2, 3,4. 

J-jj 

(2.18) 

Then we can get the following analyticity and boundedness properties: 


k-i ■ (T>3, 772, D 2 ), 

1^2- ( 774 , 77i, Di), 

hs • ( 772 , 773 , 1 ) 3 ), 

fif : {Di, D4, D4). 


(2.19) 
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2.6. Symmetries. We will show that the eigenfunctions /ij(T, t, k) sat¬ 
isfy an important symmetry. 


Lemma 2.2. The eigenfunction t, k) of the Lax pair (El satisfies 
the following symmetry: 

= A^{x^t,k) A, (2.20) 


where 

/ 1 0 0 \ 

A= 0 -a 0 , = 1. 

\ 0 0 -a J 

Here, the superscript T denotes a matrix transpose. 


( 2 . 21 ) 


Proof. The equation (12.201) follows from the fact 

— AU{x,t,k)A= U{x,t,k)'^, —AV{x,t,k)A = V{x,t,kY', (2.22) 


and 


—U{x,t,k)'^'^^{x,t,k), '^f{x,t,k) 


—V{x, t, k)'^'^^{x, t, k) 
(2.23) 

□ 


Remark 2.3. This lemma implies that the eigenfunctions pLj{x,t, k) of 
Lax pair H2.^) satisfy the same symmetry. 

2.7. The Jm,n^ computation. Let us dehne the 3 x 3—matrix value 
spectral functions s{k), S{k) and pL^k) by 

/i3(x, t, k) = p, 2 {x, t, (2.24a) 

pi{x, t, k) = /i2(T, t, ^2.24b) 

pii{x, t, k) = ixsix, t, ^2.24c) 


s(fc) = p3(0,0,fc), (2.25a) 

S{k) = pi(0, 0, k) = T, k), (2.25b) 

S^ik) = pA{L,0,k) = e-^^’^"^^pf\L,T,k). (2.25c) 


Thus, 
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And we deduce from the properties of pj and fif that S{k), 

and {s^{k), S^{k), S^{k)} have the following boundedness properties: 


s{k) : 

i^D^ U 7^4, Di U D 2 , Di U D 2 

S{k)-. 

(D 2 U D 4 , Di U Di U 

Siik) : 

{D 2 U D 4 , Di U D^, Di U 

s^(A:) : 

{Di U D 2 , D 3 U D 4 , D 3 U D 4 

S^{k) : 

{Di U D^, D 2 U 744 , D 2 U 7^4 

Sm : 

(77i U 7^3, D 2 U 774 , D 2 U D 4 

Moreover, noticing that 



Mr,{x,t,k) = k e Dr,. (2.26) 


Proposition 2.4. The Sn can be expressed in terms of the entries of 
s{k),S{k) and Siik) as follows: 


Si = 


/ 

1 

-Ai 2 

Ai 3 

\ 

/ 5 ii 

SI2 

Sl 3 

\ 

mil (A) 





0 

-A22 

A23 


, ^2 = 

S21 

S22 

■S23 




(STsA)rr 


V 

0 

A32 

•^33 

/ 


Ssi 

S32 

■S33 



(2.27a) 




^4 = 



m 33 {s)m 2 i {S)-m23 {s)m 3 i (S) 


( 5 -^ S^)ii 

77433 (5)mii (5)-mi3 (s)m3i (S') 

( 5 -^ S^)ii 

77432 (s)mii(S)-mi 2 ( 5 ) 77431 ( 5 ) 

m23 (s)mii (5)-mi3 (s)m 2 i (S) 

(s^ 

m 22 {s)mii {S)-mi 2 {s)m 2 i (S) 

0 

m33{A) 

All 

{s^'S^)ll 

0 \ 
m32 (A) 

All 


WS^l / 

m23(A) 

All 

m22 (A) , 

All / 


(2. 


where A = is a 3 x 3 matrix, which is defined as A = 

s{k)e~^^^^SL{k). And the functions 


{S'^s^)ii = S'iimii(s) - 5'2im2i(s) + S'3im3i(s), 
(s'^S'^)ii = siimii(S') - S2im2i(S') + S3im3i(S'). 


Proof. Firstly, we dehne Rn{k),Tn{k) and Qn{k) as follows: 

Rr,{k) = e-^^^"^'^Mn{^,T,k), (2.28a) 
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Qn{k) = k). 

Then, we have the following relations: 

Mn{x,t,k) = fii{x,t, k)e^^^^~^‘^'^’^^^^^Rn{k), 

Mn{x,t,k) = fi2{x,t, 

Mn{x,t,k) = /i3(a;, t, 

Mn{x, t, k) = /i3(a;, t, 

The relations fl2.29p imply that 

s{k) = Sn{k)T-\k), 

S{k) = SMR-\k), 

A{k) = S^{k)Q-\k). 

These equations constitute a matrix factorization problem which, given 
{s(fc), S'(fc), Siik)} can be solved for the {Rn, Sn, Tn, Qn}- Indeed, the 
integral equations (12.91) together with the dehnitions oi {Rn, Sn,Tn,Qn} 
imply that 

iRn{k))ij = 0 if 7"- = 7i, 

{Sn{k))ij = 0 if 7"- = 72, 

{Tn{k))ij = 6ij if 7". = 73, 

{Qn{k))ij = 6ij if 7". = 74. 

It follows that fl2.30p are 27 scalar equations for 27 unknowns. By 
computing the explicit solution of this algebraic system, we arrive at 

dZSip. □ 

Remark 2.5. Due to our symmetry, see Lemma \2.A obtained in the 
above subsection we can replace the minors by conjugate terms among 
the representation of the functions Sn{k). It may looks like much simple 
to compute the jump matrices Jm,n{x,t, k). 

2.8. The residue conditions. Since/i 2 is an entire function, it follows 
from fl2.26p that M can only have sigularities at the points where the 
S'^s have singularities. We denote the possible zeros by {kj}^ and 
assume they satisfy the following assumption. 


13 

(2.28b) 

(2.28c) 


(2.29) 


(2.30) 
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Assumption 2.6. We assume that 

• mii{A){k) has uq possible simple zeros in Di denoted by {kj}^°; 

[k) has ni — Hq possible simple zeros in D 2 denoted by 

{ kjVnUi ; 

• (s^S'^)ii(fc) has 712 —17-1 possible simple zeros in D 3 denoted by 

{kjVnUll 

• Aii{k) has N — n 2 possible simple zeros in denoted by 

and that none of these zeros coincide. Moreover, we assume that none 
of these functions have zeros on the boundaries of the ’s. 


We determine the residue conditions at these zeros in the following: 


Proposition 2.7. Let {Mn}\ be the eigenfunctions defined by H2.9\) and 
assume that the set {kj}^ of singularities are as the above assumption. 
Then the following residue conditions hold: 


Resk=kj[M]i 


A.33(fej)[M(fcj)]2 - A22,{kj)[M{kj)\2, ^2e{ki) 
mii{A){kj)m 2 i{A){kj) 


Resk=kj 


[M] 


_ S2i(kj)s33{kj)-S3i(kj)s23ikj) ^ 2 e(ki) \1 

{S^sA)33(k,)m,,{kp ^ 

I S3l{kj)s22{kj)-S2l{kj)s32(kj) „2e(ki)l ]\j:u 11 
{S^s^)33{kj)mMkp ^ 

no + l < j < ni, kj e D 2 , 


1 < i < no, kj e Di 
(2.32a) 

(2.32b) 


Resk=kj[M ]2 


m33{s){kj)M2i{S)ikj)-m23{s)ikj)M3i{S)(kj) ^- 2 B{kP\-h/r t U M 

- AsAaAAA) - 

ni + l<j< n 2 , kj E D^, 

(2.32c) 


Res - \M] = {s)(kj)M 2 i{S){kj )-m 22 {s){kj)M3i{S)(kj) 2 e{kP r ^ 11 

{sTsAwAiiikj) ^ ^ RR 

ni + l < j < 772 , kj E D 3 . 

(2.32d) 

ReSk=k,[M]^ = n 2 +l < j < N,kj E D,. 

^11 \^j) ^21 \^j ) 

(2.32e) 
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ReSk=k,[M]s = n^+l < j < N,kj e D,. 

■^11 \^j) ^21 \^j ) 

(2.32f) 

where f = and 9 is defined by 


9{x, t, k) = ikx + 2ik‘^t. 


(2.33) 


Proof. We will prove fl2.32a|) . fl2.32cp . the other conditions follow by 
similar arguments. Equation fl2.26p implies the relation 

Ml = (2.34a) 


Ms = (2.34b) 

In view of the expressions for Si and S's given in fl2.27p . the three 
columns of fl2.34ap read: 

1 


[Mi]i - [p 2 ]i- , 

mii[A) 

[Mi]2 = [ja-fiie Ai2 + [/i2]2-4,22 + [p2]3‘^32, 
[^1)3 = [h2]ie ^^^13 + [Ai2]2-4.23 + [At2]3‘^33. 
while the three columns of fl2.34bp read: 

[Mfi\i = [p2]lSll + [p2]2S2ie^^ + [p2]3-S3ie^^ 
[M3]2 = [^^2] 1 

I r,. 1 m33(s)mii(S)-mi3(s)m3i(5) 

+ [/l 2 j 2 - 


^2 3 - 


_n 

'W^ii 


I \k-z »^ 23 (^)"^ll( 5 ')-rrai 3 (s)m 2 i( 5 ) 


[M 


3J3 


rn^2{s)rn2i{S)-rn22{s)rnz\{S) ^-20 

11 


- [h2] 1 ■ ■ 11' e 

+ [h2j 2- {s^'sAll - 

_j_ __ rn,22 ( 3 )mii (S)-mi2 {s)m2i (S) 


(2.35a) 

(2.35b) 

(2.35c) 

(2.36a) 

(2.36b) 

(2.36c) 


^2J3 

We hrst suppose that kj G Di is a simple zero of mii(^)(A:). Solving 
fl2.35bp and fl2.35cp for [p 2 ]i, [^ 2)3 and substituting the result in to 
fl2.35al) . we hnd 


_ >^33 [^ 1)2 - > 4.32 [Ml ]3 

^ mii(^)m 2 i(^) 


[^ 2)2 ^2e 


m2i{A) 
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Taking the residue of this equation at kj, we hud the condition (I2.32ap 
in the case when kj G Di. 

In order to prove fl2.32cp . we solve fl2.36ap for [n 2 ]i, then substituting 
the result into fl2.36bD and fl2.36cp . we hnd 


[Msh 

[Mals 


m33[s) 

Sll 


m32{s) 

Sll 


[h: 


'2 2 - 


i«23(s)r 1 m33{s)m2i{S) - m23{s)m3i{S) _2ein^ ) 


[h: 


Sll 


^22(5 


2 2 - 


Sll 




2 J 3 - 


(2.37a) 

^32(5)77121(5') - m22(s)m3i(S') 20 


(s^S'"^)iiSii 


[M 


3jl 


(2.37b) 

Taking the residue of this equation at kj, we hnd the condition fl2.32cp 
in the case when kj E D 3 . □ 


2.9. The global relation. The spectral functions S{k),SL{k) and 
s{k) are not independent but satisfy an important relation. Indeed, 
it follows from fl2.24p that 

t, k)e^^^^+^^’‘'^^^{S-\k)s{k)e-^’‘^^SL{k)} = t, k). (2.38) 

Since /ii(0,T, fc) = I, evaluation at (0,T) yields the following global 
relation: 

S-^{k)s{k)e-^^^~^SL{k) = e-2*^'^^c(T, k), (2.39) 

where c{T, k) = /i4(0, T, k). 


3. The Riemann-Hilbert problem 

The sectionally analytic function M(x, t, k) dehned in section 2 sat- 
ishes a Riemann-Hilbert problem which can be formulated in terms of 
the initial and boundary values of qi{x,t) and q 2 {x,t). By solving this 
Riemann-Hilbert problem, the solution of fll.ip can be recovered for all 
values of x, t. 

Theorem 3.1. Suppose that qi{x, t) and q 2 {x, t) are a pair of solutions 
of / li.il) in the interval domain R. Then qi{x,t) and q 2 {x,t) can be re¬ 
constructed from the initial value {qio{x), q 2 o{x)} and boundary values 
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{ 9 oi{i), 9 o 2 {t),gii{t),gi 2 {t)}, {f 01 (t), / 02 (^), /ii(i), / 12 {t)} defined as fol¬ 
lows, 


qw{x) = qi{x,t = 0 ), q 2 o{x) = q 2 {x,t = 0 ), 

9 oiit) = qfix = 0,t), go 2 {t) = q 2 {x = 0,t), 

foi{t) = qi{x = L,t), = q 2 {x = L,t), (3.1) 

9ii{t) = qix{x = 0,t), gi2{t) = q2x{x = 0,t), 

= qix{x = L,t), /i 2 (t) = q2x{x = L,t). 

Use the initial and boundary data to define the jump matrices Jm,n{x, t, k) 
in terms of the spectral functions s{k) and S{k), SL^k) by equation 

Assume that the possible zeros {kj}^ of the functions mii{A){k), 
{S'^s^)n{k), S^)u{k) and Aifik) are as in assumption \2.6[ 

Then the solution {qi{x,t),q 2 {x,t)} is given by 

qi{x, t) = 2i lira {kM{x, t, k))i 2 , <? 2 (t, t) = 2i lim {kM{x, t, k))i 3 . 

k^oo k^oo 

(3.2) 

where M{x, t, k) satisfies the following 3x3 matrix Riemann-Hilbert 
problem: 

• M is sectionally meromorphic on the Riemann k—sphere with 
jumps across the contours DnUDm-, n,m = 1, • • • ,4, see Figure 

m 

• Across the contours Dn fl D^, M satisfies the jump condition 

Adn^x,t,kjj Adyyii^x, t, k^ JrfYi^Yii^x, t, kfi k G D,ri n D^n, n, m 1, 2, 3,4. 

(3.3) 

• M{x,t,k) =1 +Ofidfi k ^ 00 . 

• The residue condition of M is showed in Proposition |il. T\ 


Proof. It only remains to prove fl3.2|) and this equation follows from 
the large k asymptotics of the eigenfunctions. □ 
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4. Non-linearizable Boundary Conditions 

A major difficulty of initial-boundary value problems is that some of 
the boundary values are unkown for a well-posed problem. All bound¬ 
ary values are needed for the dehnition of S{k), SL^k), and hence for 
the formulation of the Riemann-Hilbert problem. Our main result ex¬ 
presses the unknown boundary data in terms of the prescribed bound¬ 
ary data and the initial data via the solution of a system of nonlinear 
integral equations. 

4.1. Asymptotics. An analysis of fl 2 . 2 p shows that the eigenfunctions 
{Hj}\ have the following asymptotics as A: —)■ cx): 


^j{x,t,k) = 1 + ^ 


■fc 2 


/ fOO A 




= 1 + 4 




kl2 

+S3 

\ 


/ 

(2) 

+11 

(2) 

+12 

(2) 

+13 

\ 

,,(b 

,,(b 


1 1 


(2) 

,,(2) 

,,(2) 


+22 

+23 


fc2 


+21 

+22 

+23 


+32 

+S 

) 


V 

(2) 

+31 

(2) 

+32 

(2) 

+33 

/ 

gi 

2i 



92 

2i 


\ 





+ 


/ 


V 


crqi 

It1^^22 

r(x,t) 

2i 

l(xj,tj) 

(yq2 


r(x,t) 

2i 

l(xj,tj) 

(2) 

+11 

(1) , (1) 
91/^22 +'12^32 

2i 

+ jQh 


A(i) 

^23 

A(i) 
^33 


/ 


V 


- 


— n 

2iQlkll 

— n 

2 i?2hil 


( 2 ) 

k22 

( 2 ) 

h32 


( 1 ) , ( 1 ) 
gi/^23 

2i 

( 2 ) 

+23 

( 2 ) 

+33 


+ \(l23 


\ 


/ 


+ 0{^). 
(4.1) 




An = cr[|(|gip + \q2\‘^)dx + |(gign - QiQix + q2q2x - q2q2x)dt\ 
aS = CT[-^\qi\‘^dx - liqixqi - qiqix)dt] 

A 23 = cr[-^qiq2dx - l{qixq2 - qiq2x)dt] 
aS = (r[-^q2qidx - \{q2xqi - q2qix)dt] 

AS = (7[-'2\Q2?dx - \{q2xq2 - q2q2x)dt]. 

( 4 - 

The functions = 1,2 are independent of k. 
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Remark 4.1. Became we do not need the asymptotic expressions of 
the and {t^if}ij =2 following analysis, we do not write down 

the explicit formulas for these functions. 


We define functions {^ij(t, and {(fijit, by: 


/i 2 ( 0 ,t, fc) = 


^ ^u{t,k) ^i2it,k) 4)i3(t, fc) ^ 

^2i{t,k) ^22{t,k) ^23{t,k) 

\^3l{t,k) ^32{t,k) ^33{t,k) J 


^ (j)u{t,k) (j)i2it,k) 013(f, A:) ^ 

ri3{L,t,k)= (j)2l{t,k) ((>22{t,k) 023(0^) 

\03l(t, A:) (j)32{t,k) (j)33{t,k) J 

From the asymptotic of fij{x,t, k) in fl4.1l) we have 


/X2(O,0/c) = 



/ 


(t) 


it) 

ifs’ 

it) 

\ 

+ 

1=1 


4<i' 

(t) 


it) 

■!>£’ 

it) 



V 

^31 

(t) 

^32 

it) 

^33 

it) 



/ 

Si'll 

it) 

$( 2 ) 

^12 

it) 

$( 2 ) 

^13 

it) 

\ 

+ F 


^21 

it) 

$( 2 ) 

^22 

it) 

$( 2 ) 

^23 

it) 



V 

^31 

it) 

<T)*^ 2 ) 

^32 

it) 

<T)( 2 ) 

^33 

it) 

/ 


+ 


(4.3) 


(4.4) 


(4.5) 


Recalling that the dehnition of the boundary data at a; = 0, we have 

^ -^12 (t) + h)^23^ (b+902 (b) 


^u{t) = 

(^) = ^f^02(t), 

®S(^) = -f/o(^ii(^)^oi(^) - goiit)gii{t))dt, 
^32 (t) = -f /o(^ 12 (t)^Ol(^) - ^ 02 (t)^ll(t))dt, 


•^23 (^) = -f /o(^ll(^)^ 02 (^) - g 0 l{t)gi 2 {t))dt, 

= -f /o(^ 12 (t)^ 02 (t) -^ 02 (t)£/i 2 (t))dh 
(4.6) 


In particular, we hnd the following expressions for the boudary values 
at X = 0 : 


^01 (t) = 2i^^.]{t), go 2 {t) = 2i^f^{t), (4.7a) 


5'ii(^) — 44*12^(t) + 2z(5(oi(t)4>22^(t) + go2{t)^f2 (t)) ^ 
gi2{t) = 44>S5(t) + 2z(^oi(f)4>S(^) + ^02(t)4>^^ (t)) 


(4.7b) 
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Similarly, we have the asymptotic formulas for t, k) = 


fi3{L,t,k) = 



/ 

A 

(t) 

0 S 2 

it) 

‘/'Ss 

it) 

\ 

i + i 


^l’ 

it) 


it) 


it) 



V 

'(’m 

it) 


it) 


it) 



/ 

dt 

it) 


it) 

4? 

it) 

\ 

+ F 


4? 

it) 

0® 

it) 

4? 

it) 



V 

4? 

it) 

4? 

it) 


it) 

/ 


(4.8) 


+ 0(4)- 


Recalling that the dehnition of the boundary data a.t x = L, we have 

= i/o2(i), 4?(i) = + Ife. W4'10^.2 (<*w) , 

- foi{t)fu(t))dt, (ji'ilit) = -f/„‘(/ii(*)/i)2W - foiit)fi2(t))dt, 


nJW = -f/„(/i2(i)/oi(0 - fo 2 {t)fn{t))dt, 4>y{t) = -f/„(/i2(i)/o2(«) - fta(t)fi 2 {t))dt. 

(4.9) 

In particular, we hnd the following expressions for the boudary values 
at X = L: 

/oi (t) = {t ), /o 2 (t) = (t ), 

fn{t) = ^uit) + + /O2(t)032 W)) 

fuit) = 40® (f) + 2i(/oi(t)053^(f) + fo2it)(p33 (t)) 

From the global relation (12.391) and replacing T by t, we hnd 


(4.10a) 

(4.10b) 


P 2 ( 0 , t, k)e^^’^"^^{s{k)e-^^^^SL{k)} = c{t, k) 


—ikLA ( 


(4.11) 


From the relation fl2.25cl) and the symmetry fl2.20p . we know that 
the spectral function 5 'l(/c) can be expressed by {0p(t, So if 

we denote the matrix-value function c{t,k) as c{t,k) = {cij{t, k))^ 

The functions {cij{t, k)}^^ij ^2 analytic and bounded in Di away 
from the possible zeros of mii{A){k:) and of order ask ^ oo. 

In the vanishing initial value case, the asymptotic of cij{t, k),j = 2, 3 
becomes much more simple. 


Lemma 4.2. We assuming that the initial value and boundary value 
are compatible atx = 0 andx = L (i.e. atx = 0, gio(O) = 5'oi(0), ?2o(0) = 
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^ 02 ( 0 ); at X = L, qio{L) = foi{0),q2o{L) = /o2(0)j. Then, in the van¬ 
ishing initial value case, the global relation ^-11^ implies that the large 
k behavior of Cij(t,k),j = 2,3 satisfies 


Ci2{t, k) 


cifit,k) 


k ' k‘^ 

I I nf 1 I 

k ' fc 2 ' ) 


+ ^(f) 

g 2 ifcL ^ 


(4.12a) 


_ I ^13 (*)+^12 (*)42 (*)+'^13^ (*)^.33 (*) 

k ' 

^31^ M I ^fi h) I nf 1 1 

k ' IP ' '^IP ' 


+ ^(f) 

g 2 ifcL /j. 


(4.12b) 


Proof. The global relation shows that under the assumption of vanish¬ 
ing initial value 

Cuit, k) = <l>i2(t, k)(j)22{t, fc) + ^isit, k)(j)23{t, k) - a^ii(l)2i{t, 

(4.13a) 

Cifit, k) = <hi2(t, fc)032(t, ^ + *hi3(t, k)fi^fit, k) - a^u(j)3i{t, 

(4.13b) 

Recalling the equation 

fit - 2ik‘^[A, fi] = V2fi. (4.14) 

From the first column of the equation (I4.14p we get 

{ 4*114 = 2A:(5'oi4>2i -|- fi'o24’3i) — *c’‘(|fi'oiP + |5'02p)4*ii -|- f(5'ii4>2i -|- 5'i24*3i), 

4*214 — 4iA:^<l*2i = 2crA:^oi4*ii -|- *c’'(|g'oip4>2i + ^ 015 * 024 * 31 ) — zcr^ii4*ii, 

4*314 - 4iP4>3i = 2akgo2^u + icr(^o25*oi4*2i + |5o2p4>3i) - iagi2^n, 

(4.15a) 

From the second column of the equation (14.141) we get 

{ 4*124 + 4iA:^4>i2 = 2fc(5foi4*22 + 5'024*32) — ZCxd^'oiP -|- |5'02p)4*l2 + «(5ll4>22 + 5'124*32), 
4*224 = 2 (T/c^oi 4*12 + *<T(|5*01p4>22 + ^0l5024*32) “ *(L5ll4>i2, 

4*324 = 2akgo2^12 + *<t(^0250i 4*22 + |5*02p4>32) — iO'gi2^12, 


(4.15b) 
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From the third column of the equation (14.141) we get 

*^134 + 4iA:^$i3 = 2k{goi^23 + 9'02‘h33) — ^u’d^'oiP + |5'02p)‘hi3 + i(5'll*h23 + 5'l2‘h33), 
*^23* = 2(Tfc^oi‘hi3 + io'{\goi\‘^^23 + ^oi5'02‘h33) — iagn^i3, 

*^*33* = 2(TA;^02*hi3 + *cr(^025'0l‘^*23 + Ifi'02p*h33) — Z(T()i2*hl3, 

(4.15c) 

Suppose 


*^21 

\ *^31 / 






+ 


fc2 


+ • • ■ e 


(4.16) 

where the coefficients aj{t) and j = 0,1, 2, • • •, are independent 

of k and are 3x1 matrix functions. 

To determine these coefficients,we substitute the above equation into 
equation fl4. 15al] and use the initial conditions 

Q^o(O) + /3o(0) = ( 1 0 0 )'^, cti(O) + /3i(0) = ( 0 0 0 )'^. 


Then we get 

^’21 

V <^31 / 


/ 1 \ 
0 

Vo/ 




<F 


11 

( 1 ) 

21 




+ 


fc2 


4>ii’ / 




0 \ 


+ 

1 

k 

-bi'(o) 

+ 0{^) 



-bi'(o)) 





<F 


11 

( 2 ) 

21 


■S>S’ / 


0(i 




(4.17) 


Similarly, suppose 
/ $12 \ 


*^22 
V *^32 / 


= (a„(i) + 2:^ + ^ + ...)+(^A(() 




k 

(4.18) 

where the coefficients aj{t) and j = 0,1, 2, • • •, are independent 

of k and are 3x1 matrix functions. 
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To determine these coefficients,we substitute the above equation into 
equation fl4.15bp and use the initial conditions 


Q^o(O) + /5o(0) — ( 0 1 0 )^, tti(O) + /5i(0) — ( 0 0 0 


Then we get 


/ <^>12 \ 
*^*22 

\ *^*32 / 


/0\ 




«!? \ 

1 

+ 1 


+ F 

4g> 



) 


■»!>? / 


Oil 






/ 

-kS 

(0) \ 


/ 

-$(5(0) + $S5(o)<^>S + <^>S(o)<f>S 

\ 


+ 

1 

k 


0 


+ F 


-f^oi$J5(0) 


+ o(^) 



V 

0 

) 


V 

-f^02$S5(0) 






(4.19) 

Similar to the derivation of ^i 2 ,i = 1,2,3, from fl4.15cl) we can get 
the asymptotic formulas of 4>i3, z = 1, 2, 3 



/0\ 

0 

V 1 / 


+ \ 


(\ 
V / 






+ 0{h) 




/ 

-kS 

(0) \ 


( 

-KfhO) + 1>Sho)<l>S + -S'lhojl’ss 

\ 


+ 

1 

k 


0 


+ F 


-fSo.tSUo) 


+ 0{^) 



V 

0 

) 



-f9o2>I>Sh0) 




^—Aik^t 


(4.20) 

Similar to (14.151) . we have the equations {(l>ij}%j=i satisfy the similar 
partial derivative equations: 

From the hrst column of the equation fl4.14p we get 


hit = 2 A:(/oi 021 + /o2</>3l) ~ *Cr(|/oiP + |/o2p)011 + *(/ll^21 + /l2<(>3l), 
hit — 4z/c^02i = 2(tA:/oi0ii + zcr(|/oip02i + foifo2hi) ~ 

hit - = 2akfo2hi + *cr(/o2/oi02i + I/02I Vsi) - icrfuhi , 

(4.21a) 
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From the second column of the equation (14.141) we get 

{ 4‘l2t + 4ifc^0i2 = 2 A:(/oi 022 + /o2032) ~ *Cr(|/oiP + |/o2p)012 + *(/ll^22 + /l2</>32) 
022 i = ‘^Crkfoi(j)l2 + *0'(|/oiP^22 + /oi/o2032) — ^Cr/ll012, 

032i = 2 ctA:/o 2012 + *0'(/o2/oi<?^22 + |/o2p032) ~ *<7/12012; 

(4.21b) 

From the third column of the equation fl4.14l) we get 

{ 013i + 4ifc^0i3 = 2 A:(/oi 023 + /o2033) “ *<7(|/oiP + |/o2p)013 + *(/ll023 + /l2033) 
023i = 2crA:/oi0i3 + *<7(|/oiP023 + /oi/o2033) “ *<7/ll013; 

0331 = 2cr/i:/o2013 + *<7 (/o2/o 1023 + |/o2p033) “ *<7/l2013; 

(4.21c) 

Then, substituting these formulas into the equation fl4.13ap and 
noticing that we assume that the initial value and boundary value 
are compatible at x = 0 and x = L, we get the asymptotic behavior 
04.12 all of Cij{t, k) as k ^ oo. Similar to prove the formula 04.12bp . 

□ 

4.2. The Dirichlet and Neumann problems. We can now derive 
effective characterizations of spectral function S{k), S^lk) for the Dirich¬ 
let ({5'oi(^),5'02(t)} and {/oi(f),/o2(t)} prescribed), the Neumann ({t/ii(t), t/i2(t)} 
and {/ii(t),/i2(^)} prescribed) problems. 

Dehne functions as 


f-{t,k) = f{t,k) - f{t,-k), f+{t,k) = f{t,k) + f{t,-k), (4.22) 

Introducing 

A{k) = S(A:) = (4.23) 

Denoting dD^ as the boundary contour which is not included the zeros 
of A(fc). 

Theorem 4.3. LetT < oo. Let qio{x),q2o{x),0 <x <L, he two initial 
functions. 

For the Dirichlet problem it is assumed that the function {goiit), go 2 it)}, 0 < 
t < T, has sufficient smoothness and is compatible with qio{x), q 2 o{x) 
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at X = t = 0, that is 


9io(0) — 5'oi(0), g 2 o( 0 ) — g'o 2 ( 0 ). 

The function {foi(t), fo2(t)},0 <t <T, has sufficient smoothness and 
is compatible with qio{x), q 2 o{x) at x = L, that is, 

qio{L) = /oi(0), q 2 o{L) = /o 2 ( 0 ). 


For the Neumann problem it is assumed that the function gi{t),0 < 
t < T, has sufficient smoothness and is compatible with qo{x) at x = 
t = 0; the function fi(t),0 < t < T, has sufficient smoothness and is 
compatible with qo(x) at x = L. 

For simplicity, we suppose that mii{A){k) has no zeros in Di. 

Then the spectral function S{k) is given by 


$11 (A:) \ 

S{k) = I 


-a$i3(/c)e 


-Aik^T 


^22{k) 

^23{k) 


^32{k) 

$33 (fc) J 


( 


Srik) = 


Mk) 

-a(j)i2{k)e~‘^^’^^'^ 

-Aik^T 


(4.24) 


222 


Ck) 


\ -a(f)i 3 {k)e 


(p23{k) 


<f32{k) 

f)33(k) 


! 

(4.25) 

and the complex-value functions {^is^t, /c)}f=i satisfy the following sys¬ 
tem of integral equations: 


^i3{t,k) = /o e + ko2n$i3 

+ {2kgoi + ign)^23 + (2%o2 + *5'i2)$33] k)dt' 


(4.26a) 


$23(A, k) — a [{2kgoi — ign)^i 3 + i\goi\‘^^23 + *5'oi5'02$33] {t', k)dt' 

Jo 

(4.26b) 

$33(A) k) = 1+ [ a [( 2 A :^02 — *^12)$13 + *5'01^02$23 + d5'02|^$33] (A^ k)dt' 

Jo 

(4.26c) 
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and satisfy the following system of integral 

eguations: 


k) — 1 + Jp [—^o'dS'oiP + |fl'02p)‘^’ii 

+ {2kgQi + igii)^2i + ( 2^02 + *5'i2)*^’3i] k)dt' 


(4.27a) 


$21 (t, k) — [ * V [(2A;5foi ~ + dfl'oiP ‘^21 + *5'oifl'o2*^’3i] (^^ k)dt' 

Jo 

(4.27b) 

*^’3i(d k) = [ * V [(2A;^02 — *^i2)‘^’ii + *fi'oi^o2‘^’2i + d5'02|^‘^’3i] (^^ k)dt' 

Jo 

(4.27c) 


^i2{t,k) =/oe + Ifl'02n<^>i2 

+ (2^01 + *5'ii)‘^22 + (2%o2 + ^5'i 2)‘^’32] (^d k)dt' 


(4.28a) 


$22(d ^) — 1+ / C’’ [(2 /i:5'01 ~ *5'1i)*^’ 12 + d5'0lp‘^22 + *fi'0lfi'02*^’32] {t\ k)dt' 
Jo 

(4.28b) 

‘^ 32(^5 k) = [ O' [{2kgo2 — igi2)^12 + *fi'0l'702‘^’22 + dfi'02p*^’32] {t', k)dt' 

Jo 

(4.28c) 

Funetions { 4 >ij{t,k)}F^^ satisfy the same integral equations replaying 
thefunetions {goi, go2, gii, gu} with {/oi,/ 02 ,/ii,/i 2 }. 

(i) For the Dirichlet problem, the unknown Neumann boundary 
value {gu{t),gi2{t)} and {fn{t)Ji2{t)} are given by 

^ii(0 = ^ IdDi f (^^12- + igoi)dk + f /g^o(^oi<^>22- + go2^32-)dk 

fdDo(9oi(^22- + go 2 (j) 23 -)dk - ^ Jg^o ^{k(f) 2 i- + iofoi)dk 

fano Z [(*^12(022 — 1) + ‘^*13023) + (j($ii — l)(^2i] _ dk. 

(4.29a) 

gi2{t) = ^ fgjjo §(A:$ 13 - + ig02)dk + | /g^o(fl'Ol‘^’ 23 - + g02^33-)dk 
~n IdD°^ 9 oi 4 >^ 2 - + go 2 <p 33 -)dk — ^ Jgjgo -^(/s^si- + iof02)dk 

fdD° Z [(‘^12<^32 + 4>i3((^33 — 1 )) + (j($ii — l)(^3l] _ dk. 

(4.29b) 
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and 

/ll(^) = —^ + ifoi)dk + ^ /g^o(/oi022- + fo24>32-)dk 

fdDl(foi^22- + fo2^23-)dk + fgj^o ^(k^2i- + 'icrgoi)dk 

IdD^ Z “ 1)5*21 — (012(^22 ~ 1) + 0135*23) _ dk. 

(4.30a) 

fuit) = —^ fgjjo ^(^ 013 - + ifo2)dk + ^ /g^o(/o 3023 - + /o2 033 -)c^^ 
fdDl(foi^32- + fo2^33-)dk + ^ fg^o ^(k^3i- + icrgo2)dk 

fdD° Z [*^(011 “ 1)5*31 — (0125*32 + 013(5*33 “ 1)) _ dk. 

(4.30b) 

where the conjugate of a function h denotes h = h(k). 

(ii) For the Neumann problem, the unknown boundary values {goi(t), go 2 (t)} 
and {foi{t), fo 2 {t)} are given by 

9oi{t) = i jgjg, ^^i2+dk - I fg^o ^^21+dk 

~w IdD^ K ['^(‘1*11 “ 1)021 ~ (^12 (022 ~ 1) + ^13023)e ^ dk, 

(4.31a) 

902{t) = i fg^o i^i3+dk - ^ fg^o z^3i+dk 

~w IdD^ 'k ['^(‘1*11 “ 1)031 ~ (*1*12 032 + *1*13 (033 ~ 1))^ ^®^^] ^ dk, 

(4.31b) 

and 

ki{t) = -^ fgjjo if>i2+dk - § fg^o ^^21+dk 

fdD° k [^(011 “ 1)5*21 — (012(^22 ~ 1) + 0135*23)e^*^^] ^ dk, 

(4.32a) 

fo2{t) = fg^o ik3+dk - I fg^o ^^31+dk 

~wIdD0 k [^(011 “ 1)5*31 — (012^32 +013(^33 “ l))e^®^'^] ^ d/c, 

(4.32b) 

Proof. The representations fl4.24p and fl4.25p follow from the relation 
S{k) = e“^®^^‘^/i^^(0, t, A;) and S'l(A:) = e“^®^^*^/ij^(0, f,/c), respec¬ 
tively. And the system (??) is the direct result of the Volteral integral 
equations of /i 2 ( 0 , t, k) and jaz^L, t, k). 
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(i) In order to derive (I4.29ap we note that equation (I4.7bp ex¬ 
presses gn{t) in terms of <I) 22 ^(t, A:), <h 32 ^(t, fc), A;). Fur¬ 

thermore, equation fl4.5p and Cauchy theorem imply 


k) - l]dk = /g^J<F2x2(t, k) - l]dk, 
W = IdD^ ^32(C k)dk = ^32(t, k)dk, 

and 



k^u{t,k) 


9oi{t) 

2i 


dk 



k^i2{t,k) 


9oi{t) 

2i 


dk, 


Thus, 

i7r4>g(t) = - + fg^J [<I>22(t, k) - l]dk 

= (fdDi + IdDs) [*^22(C k) - l]dk 

= fgj:>J^22(t, k) - l]dk - /g^3[<h22(C -k) - l]dk 
= fdD3(^22(t, k) - ^22(t, -k))dk (4.33) 

= jgjg^^22-{t,k)dk 

i7r<Fg(t) = - + fg^J [4)32(t, k)]dk 

= IdDs "^32-(t,k)dk 


Similarly, 




\IdD 3 4 “ IdDi 
= IdDa 

= IdDi { 


k<l> 


k<^>u{t,k)-^-^ 

k^i2it,k) 


,2{t,k)-^-^ 

dk 


dk 


2 i 

gold) 

2 i 


+ 


2e- 


■[k^i2{t, k) 


where I{t) is dehned by 


(4.34) 



2g-2ifcL 

A 


[k^i2{t,k) 


£'oi(^) il 

2i 


dk 
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The last step involves using the global relation (??) to compute 
I{t). That is, 


m =s, 

+fi 

+ faD° 


9D° 


dDl 


2e- 


2e- 


kc^o - o - 




(1) _ ^^2't>22 't ’23 I ~jX^)p2ikL 

12 k ' Op>21 ^ 


dk 


+ cr(/c02i - 021 )e 


(l)l„2ifcL 


dk 


2e- 


■ ^k ($12(022 ~ 1) + $13 023 + Cr(‘hii — 1)0216^*^'^ 

(4.35) 

Using the asymptotic (??) and Cauchy theorem to compute 
these terms on the right-hand side of equation fl4.35p . we hnd 

I(t) = -z7r$S5 - fg^o [^022- - ^(^021- + i(jfoi)] dk 

IdD° ^ [(*^12(022 ~ 1) + $13 023) ^ 4- Cr($ii — l)02l] _ dk. 

(4.36) 

Equations fl4.33p . fl4.34p and fl4.36p together with (l4.7bH yield 
fl4.29ap . Similarly, we can prove fl4.29bp . 

The expression (Id.dOap for fu (t) can be derived in a similar 
way. Indeed, we note that equation fl4.10bl) expresses /ii(t) 
in terms of 022 ^ (t, k), 032 ^ (t, k), 012^0 These three functions 
satisfy the analog of equations fl4.33p and 14.341 In particular, 
012 ^ satisihes 

Z7r0g(t) = Jg 


IdDs 


0012 - (Pu 


dk 

(!)■ 


= IdDl I-f (0012 - </)\2 )}_dk+ J(t) 


where 


J(t) = 


'dDO 


2e 


2ikL 


A 


-(0012 - 012 ^ dk 


(4.37) 


(4.38) 


Then using the global relation (??) to compute J{t): 


m =x 


SDS 


^ 0CrC21 -f 032^6^*^-^ 




(1) I </>iy^22+<^l3^23 ^2ikL 


a$^0 + 


dk 


+ laot i \ I. * 


+ fg^O ^ [< t (011 — 1)$21 — ( 012 ( 4*22 ~ 1 ) + 013023 )^^*^^] _ dk 

(4.39) 


dk 
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The equations (14.391) . (14.371) together with the asymptotics of 
C 2 i{t, k) yield fl4.30ap . The proof of fl4.30bD is similar. 

(ii) In order to derive the representations fl4.31ap relevant for the 
Neumann problem, we note that equation fl4.7al) expresses ^'oi 
and gQ 2 in terms of and $ 13 ^ respectively. Furthermore, 
equation 04.51) and Cauchy theorem imply 

-f 4)2 (*) = iao, iuit, k)dk = 4i2((. k)dk, (4.40) 

Thus, 

= (idDs+fdDi) ^12(t,k)dk 
= fdD3^i2-(t,k)dk (4.41) 

= IdDl (§^i2+(t, k)) dk + K{t), 

where 

!<(() = - f (4.42) 

JdD° ^ 

and using the global relation and the asymptotic formulas of 
C 12 , we have 

K{t) = - fg^o |l^^2i+ + I [o'(4>ii - 1)^21 - (*hi2(022 - 1) + *hi3023e"^**'^)]_^ 

(4.43) 

Equations 04.7ap . 04.4ip and 04.431) yields 04.31ap . The proof 
of the other formulas is similar. 

□ 

4.3. Effective characterizations. Substituting into the system (??) 
the expressions 



j,0 + 

:i,2H-, i,i = 1,2,3. 

(4.44a) 

4^ij 

j,0 + ^0p',l + 

s2 + • • • , b J = Ij 2, 3. 

(4.44b) 

(1) 

^01 = £^01 

+ + • • • ) 

6^02 = ^5^02^ + ^“^9(^2 + ■ ■ ■ ) 

(4.44c) 

/oi = 

+ + • • • ) 

fo2 = + ■ ■ ■ ) 

(4.44d) 
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gii = £g[i + ^'^gii + • • •, gi 2 = ^g^ + ^'^gu + ■ ■ ■) (4.44e) 

/ll = ^fu + ^‘^fu + ■ ■ ■ 5 fl2 = £fl2^ + + • • • 5 (4.44f) 

where e > 0 is a small parameter, we find that the terms of 0(1) give 


{ *hl3,o = 0 $23,0 = 0 $33,0 = 1, 

*^11,0 = 1 $21,0 = 0 $31,0 = 0, (4.45) 

$ 12,0 = 0 $ 22,0 = 1 *^ 32,0 = 0 . 

Moreover, the terms of 0{e) give 

*^33,1 = 0 $23,1 = 0, 

$13,1 (f, k) = /g 
$11,1 = 0, 

0(e) ; < $ 21,1 = /o - ig[\^){t’)dt’ , (4.46) 

$31,1 = /o -ig\^j){t’)dt’, 

$12,1 = /o +i^ff)(t')dt', 

, $22,1 = 0, $32,1 = 0. 

and the terms of O(e^) give 

$13,2 = (2%JJ +ig^^^){t')dt’, 

$23,2 = /g a[(2A:^Ji^ -i^Ji^)(f')$i3,i(t',/i:) +i^Sl^(t')^o?(^')]c^^', 
$33,2 = /o c^[(2A:^S? +i\gl!^{t')\^]dt', 

$11,2 = Jo -i<yi\gm? + \g?‘J?)it') + J^kg^m + )(t')$2i,i(t', k) 

+ {2kgJJ+igJJ){t')^^i,i{t',k) dt', 

$21,2 = /o - i/JJ){t')dt', 

$31,2 = /o -ig^JJ){t')dt', 

$12,2 = /g + ig^JJ){t')dt', 

$22,2 = /o {‘^kg^J - ig[J){t')^i2,i{t', k) + f|^oi^(t')P dt', 

$32,2 = /o cr {2kgSJ - igJJ){t')^i2,i{t’, k) + Wm {t')gm J') dt'. 

(4.47) 

Similarly, we will have the analogue formulas for {^p,z}f,j=i, I = 0,1,2 
expressed in terms of the boundary data at x = L, that is {fij^YiZoJj ^ = 
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On the other hand, expanding (I4.29p and assuming for simplicity 
that mii{A){k) has no zeros, we hnd 

9n{t) = —f {k^i2,i-{t,k)+'i9oi)dk-— f ^{k(j) 2 i,i_+iaf^l^)dk, 

JaD° JdD° ^ 

(4.48a) 

9S{t) = —[ ik^i3,i-{t,k)+i9S2)dk-— [ ^{k(j)3i^i_+iaf^2^)dk, 

JdD° ddD° ^ 

(4.48b) 

fu\i) = -—[ ik(f)u,i-{t,k)+ifth^)dk+— [ ^{k^2i,i-+icr9Si)dk, 
7^^ Jgijo m Jq£)0 a 

(4.48c) 

/i?(^) =-— / {k(j)i3^i-{t,k)+ifo2)dk+— [ ^{k^3i^i-+icrgS^)dk, 

JanS JaD^ ^ 

(4.48d) 

We also hnd that 

*^>12,1- =4fc 
$13,1- = 4fc /o 
= 4:ak 
031,1- = 4cj/c 

The Dirichlet problem can now be solved perturbatively as follows: 
assuming for simplicity that mii(^)(A:) has no zeros and given (/om 9 'o 2 ^ 
and foWfolK we can use equation fl4.49p to determine We can 

then compute 5 'ii\g'i 2 ^ and /ii\/i 2 ^ from (??) and then $ij,i,i = 2,3 
from fl4.46p and the analogue results for 4>ji,i,j = 2,3. In the same 
way we can determine $ 1^,20 = 2,3 from fl4.47p and the analogue 
results for 0ji,2,i = 2,3, then compute 5 'ii\g'i 2 ^ and fn,fi 2 - And 
these arguments can be extended to the higher order and also can 
be extended to the systems (I4.27p . (I4.28p and (I4.26p . thus yields a 
constructive scheme for computing S{k) to all orders. The construction 
of S'i(A:) is similar. 

Similarly, these arguments also can be used to the Neumann problem. 
That is to say, in all cases, the system can be solved perturbatively to 
all orders. 
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4.4. The large L limit. In the limit L —>■ oo, the representations for 
giiii), Quit) and (7oi(^), 1702 (t) of theorem reduce to the corresponding 
representations on the half-line. Indeed, as L —)■ cxo, 


/oi —^ 0, /o2 —)■ 0, /ii —)■ 0, /i2 —)■ 0, 

4>ij ^ij, 1 as — )■ oo in Ds 


Thus, the L —)■ oo limits of the representations (I4.29ap . fl4.29bp and 
(14.3lap . fl4.31bp are 


gil{t) — ^ fgj;)o(k^l2- + Woi)dk + I /g^o(fi'Ol‘^22- + g02^32-)dk — ^ 5 'O 1022 -'^^ 

gi2{t) = ^ fgj^o(k^i3- igo2)dk | fgj^o(goi^23- + go2^33-)dk — ^ go24>33-dk. 

(4.50) 

and 


goiit) = ^ Jqdo ^i2+dk, go2{t) = ^ fg^o ^13+dk, (4.51) 

respectively, and these formulas coincide with the corresponding half- 
line formulas, see IZ21- 
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